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Abstract 

In this talk I show how to canonically quantize a massless scalar field in the 
background of a Schwarzschild black hole in Lemaitre coordinates and then present 
a simplified derivation of Hawking radiation based upon this procedure. The key 
result of quantization procedure is that the Hamiltonian of the system is explicitly 
time dependent and so problem is intrinsically non-static. From this it follows that, 
although a unitary time-development operator exists, it is not useful to talk about 
vacuum states; rather, one should focus attention on steady state phenomena such 
as the Hawking radiation. In order to clarify the approximations used to study this 
problem I begin by discussing the related problem of the massless scalar field theory 
calculated in the presence of a moving mirror. 
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1 Introduction 



The phenomenon of Hawking radiation has fascinated physicists ever since Hawking's 
1974 paper0 showed that, when quantum fields are taken into account, a black hole 
of mass M appears to emit (nearly) thermal radiation with a temperature Th = 
l/(87rGM). This result, extensively studied in the literature (see for example 0, ^, 
H), is clearly robust since all approaches lead to the same conclusion. So far as I 
know, however, no derivation of Hawking radiation discusses the problem within the 
framework of canonical quantization. This is one reason the question of whether the 
time evolution of the system is unitary became a topic of debate. 

In this talk I present work|jl|] done in collaboration with Kirill Melnikov, show- 
ing that a simple canonical quantization procedure exists and leads to the usual 
Hawking results within the framework of a unitary quantum theory. This theory 
has some surprising features but doesn't seem lead to paradoxes which invalidate 
the approach. One unique advantage of our calculation is that we can compute the 
energy-momentum tensor at all positions and times. This allows us to explicitly ex- 
hibit transient behavior and retardation effects, and - as a matter of principle - give 
a self-consistent treatment of back reaction. 

Before diving into details I will spend a few moments reminding you why carrying 
out a Hamiltonian formulation of the problem of a massless scalar field in the presence 
of a black-hole background seems so difficult. Consider the usual Schwarzschild metric 
for a black hole of mass M: 

ds^ = _(1 _ ^) de + —^2177 + r^dn^ (1) 
^ (1 - ^) 

and a massless scalar field with Lagrange density 

C = ^g[g^''d,(t>{x)d,(^{x)]. (2) 

(Note that in what follows I will set 2M = 1 to simplify the equations.) While one 
usually focuses on the singularity of the metric at r = 2M, this coordinate singularity 
is not a problem. The real issue is that in order to canonically quantize this theory we 
need a family of spacelike slices which foliate the spacetime. Given this, we can then 
define the field and its conjugate momentum on one of these slices. Inspection of the 
metric, Eq.|l], shows that surfaces of constant Schwarzschild time change from spacelike 
to timelike at the horizon (r = 2M) and so, they do not fulfill our requirements. 

2 Summary of Results 

In order to define a set of spacelike surfaces which foliate the Schwarzschild spacetime 
we begin by introducing Painleve coordinatess because surfaces of constant Painleve 
time extend from r = to r = oo and are everywhere spacelike. Since, however. 
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the Painleve timelines are not everywhere timehke we will have to make one more 
transformation, to Lemaitre coordinates, in order to canonically quantize the theory. 
In these coordinates we see that the metric and therefore the Hamiltonian is explic- 
itly time dependent, which is the reason why Hawking radiation exists. Of course, 
despite the time dependence of the Hamiltonian, the theory possesses a unitary time 
development operator, which is all we need to discuss the relevant physical issues. 
Note, because the Hamiltonian is time dependent, it is never useful to talk about 
eigenstates. 

The strategy of the computations I will describe goes as follows: first, we canoni- 
cally quantize the theory; second, we devise a scheme (the Geometric Optics Approx- 
imation) for approximately solving the Heisenberg equations of motion (this is easier 
than calculating the unitary time development operator); third, we compute the tem- 
perature seen by a thermometer located at fixed r which is adiabatically switched on 
and off; last, we compute the flux passing through a sphere of fixed radius and show 
that in the long time limit (after transients have died out) it has the Hawking value. 



3 Review of Coordinate Systems 

While we are ultimately interested in the surfaces defined by fixing Painleve time, it 
is convenient to first introduce Kruskal coordinates and plot the various surfaces of 
interest in these coordinates. The reason for this is that these coordinates make it 
particularly easy to draw nuU-geodesics (they are simply lines parallel to either the 
X OT Y axes shown in Fig.l) and they allow us to easily compare surfaces of fixed 
Schwarzschild time to surfaces of fixed Painleve time. Kruskal coordinates are defined 
by the equations: 

xy = {r-l)e' , ^\ = ' 
In these coordinates the Schwarzschild metric becomes: 

d.^ = ^^^:^+rW (4) 
r 

Note that Eq.^ immediately shows us that fixed Schwarzschild r is a hyperbola in 
the x, y-plane, as shown in Fig.l, and that surfaces of fixed Schwarzschild time cor- 
responds to straight lines x = \y\e^^ (which are not shown in Fig.l). 

Now, Painleve coordinates are arrived at by making an r-dependent shift in 
Schwarzschild time; i.e., 



t = X- 2y/¥ - In 



r — 1 



(5) 



These are the almost horizontal curves shown in Fig.l which clearly foliate the space- 
time. Note, that while Schwarzschild t and Painleve A differ by functions of r, if two 
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Figure 1: 



events having the same r then the differences in Schwarzschild and Painleve times 
are equaL In Painleve coordinates the Schwarzschild metric takes the form 

ds' = -(1 - -) dX' + + rfr^ + (6) 

Although Painleve coordinates are useful for defining our set of spacelike surfaces they 
don't work well for canonical quantization because of the cross term involving dX dr 
which appears in the metric. A better coordinate system is provided by Lemaitre 
coordinates, which are related to Painleve A and r by 

r{X,rsc.) = {rtil-lxr^' = {l{v-X)r^' (7) 

In Lemaitre coordinates the metric takes the form 

ds^ = -dX^ + —^d7]^ + r(X,r]fdn^ (8) 
r{X,r]) 

which allows a straightforward treatment of canonical quantization. 

4 Canonical Quantization 

We choose the surface A = as the initial surface on which to canonically quantize 
the free massless scalar field theory. Since the metric in all coordinate systems is rota- 
tionally invariant we can always solve the field equations for each angular momentum 
mode separately. For this reason we can imagine expanding the field in spherical har- 
monics in 6 and and then restricting attention just to the L = mode, since this 



contains most of the interesting physics. If we do this then, in Lemaitre coordinates, 
the L — scalar field Lagrangian reduces to 

>C = J kdxcl^ (A, v)? - r {d,<P (A, r,)r] (9) 



2 

where the determinant ■\/—g is 

v^ = r^/' = ^(r/-A). (10) 

Following the usual rules the momentum conjugate to the field is 

7iiX,ri) = {ri-X)d,cf>{X,ri), (11) 
and the canonical Hamiltonian is 

The commutation relations for and tt are 

[7^(A,7y),0(A,V)] = -^5(7y-7y') (13) 

As advertised, this Hamiltonian is explicitly time dependent and in a sense this fin- 
ishes our job, since we see that we are looking for steady state and not static behavior. 
Another feature of this Hamiltonian is that setting A = we see that by a simple 
rescaling of and tt, in order to absorb the factor of Ijr] in the tt^ term, we can 
convert it to the usual Hamiltonian of the L = mode of a free massless field in fiat 
space and therefore solve it exactly. 

It is important to note, as in the usual interaction representation, the fact that 
we have a time dependent Hamiltonian doesn't mean that we don't have a unitary 
time-development operator. There is a one parameter family IJ (A) which satisfies the 
equation 

^U{X) = H{X)U{X) (14) 

whose solution is the path ordered exponential of the integral of H[X). Fields at later 
A are defined by 

0(A,77) = C/(A)0(77)C/t(A) (15) 
7r(A,r^) = U{X)7r{r))U\X) (16) 

It follows from the canonical commutation relations that these operators satisfy 
Heisenberg equations of motion of the form 

dx iiv - ^)dx<P] - dr, [{v - A) r d,<P] = (17) 
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5 Solving The Heisenberg Equations: An Aside 



To explain the Geometric Optics Approximation to these Heisenberg equations it is 
convenient to first discuss these equations in fiat space. I begin with the case where 
there are no boundary conditions and then consider the problem of a moving mirror, 
whose physics is much closer to that of the black hole. By a moving mirror, I mean 
a free field theory in fiat space together with the boundary condition that the field 
vanishes on and to the left of a curve x{t). 

To solve the free field Euler-Lagrange equation when there are no boundary con- 
ditions we rewrite the equations as 

{d^-dl)<l>{t,x) = 
(dt-d,)(dt + d,)(l)(t,x) = (18) 

The general solution to this equation is 

(l)it,x)^fix-t)+gix + t) (19) 

where the functions / and g are determined by the values of 0(t, x) and its time 
derivative at i = 0; i.e., 

dtcpit = 0, xq) = -d^f{xQ) + d^g{xQ) 

a^0(i = O,xo) = d:^f{xo) + d:cg{xo) (20) 

which says that 

dxf{xo) = ^{dx(l){xo) - Trixo)) 

dxgixo) = ^{dx(l){xo) + 7r{xo)) (21) 

Now consider, as shown in Fig.2, the case of a field theory with moving boundary 
x{t) — —t + A(l — e"^***), where I have chosen to plot the curve ior A — 1 and 
K = 1/2. In this case the Euler-Lagrange equations remain unchanged, however the 
solution needs to be modified to maintain the boundary condition which says that 
4>{t,x{t)) = 0. This is easily done by the trick of adding a refiected wave, go{x — t), 
so that the general solution has the form 

(l){t,x) = e{x-t)f{x-t)+g{t + x) 

+9{t - x) go{x - t) (22) 

We will see that the crucial feature of this solution is that, if one sits at a fixed point 
X, the refiected rays contributing to (p{t,x) ior t » x all come from very near the 
point X — A. 
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Figure 2: 

6 The Story of Temperature 

In what follows I adopt Unruh's definition of a thermometer; i.e., a quantum system 
with multiple energy levels interacting with the field x). In other words, I add an 
interaction of the form 

V^i„t(t) = ee-(*-*«)'/2'^Q0(x,t) (23) 

to the Lagrangian, where the parameters to and a define the range in t for which 
the interaction is turned on and x specifies the spatial location of the thermometer. 
Furthermore, the operator Q is some operator causing transitions among the energy 
eigenstates of the thermometer. 

Of course, assumptions have to be made in order to get reasonable results. First, 
in order for the thermometer to know the mirror is moving, it is necessary that 
to >> X. Second, we impose an adiabatic condition, y/a >> 1/E, where E is the 
typical excitation energy of the thermometer, so that we do not excite the termometer 
just by turning it on or off. Finally, we impose the condition ~ k, so that the 
acceleration of the mirror is capable of exciting the higher states of the thermometer. 

With these assumptions second order perturbation theory in e tells us that the 
probability of the thermometer being excited to a state with energy E is 

V{E,Eo) = e^\{E\M{0)\Eo)\'' I dtdt' 

^-i{E-Eo){t-t')^-{t-tof/2a 

e-(*'-*o)'/2a(0(t,x)0(t',x)) (24) 

What we now do is compute {4>{t, x) 4>{t', x)) by plugging in the formula giving 0(t, x) 
in terms of (pix) and tt{x) on the original surface, expand these t = operators in 
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terms of annihilation and creation operators and evaluate the resulting expression. 
The key feature of this problem is that the points on the original surface corresponding 
to {t,x) and {t',x) for t,t' » x are exponentially close to the point x — A. A 
straightforward computation gives the result 

e'^ \{E\M{0)\Eo)\^ 
2 E — Eq 



g27r {E-Eo)/k _ I 



r{E,Eo) = 

X 



7 Computing The Energy Flux 

To compute the flux of energy through a plane at position x we need to compute Ttx, 
the energy momentum tensor for the massless field, which is given by 

^M- = ^R0(^,a;),a,0(t,a;)}, (26) 

Generally, the expectation value of components of T^i, will be divergent since the field 
derivatives are being evaluated at the same spacetime point. To deal with this we 
adopt a point splitting procedure; i.e., we define 

Ttx = I {7i{t + 6,x), dA{t - 6, x)} , (27) 

evaluate the expression and then take the limit 5^0. The result is that the energy 
density Tu diverges as but the fiux, T^^, is finite and unique. The result of the 
computation is 

Ttx = — (28) 

Although the finiteness of the flux may seem surprising, it is in fact a consequence of 
a general theorem. 



8 The Black Hole 

Now let us turn to a discussion of the case of a Schwarzschild black hole. Another way 
of describing the previous example is to say that we solve the Heisenberg equations 
by tracing back the two null-rays leaving the point {t, x) to flnd the two points at 
which they intersect the t = surface and write the field in terms of the and vr at 
those two points. There is a simple generalization of this approach to the equations 
in curved space. 

Straightforward analysis shows that what is usually referred to as a WKB analysis 
of the wave equation in curved space amounts to, what we will call, the Geometric 
Optics Approximation. The generalization is given by the following prescription. 
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First, starting from the point (A, r), find the two nuU-geodesics which meet at this 
point and trace them back to the surface A = (since this is a geometrical statement 
it can be solved in any coordinate system) . In Painleve coordinates the equations for 
these two geodesies are 

S^{x2) = X + S2{r) 
Si{r) = r-2v/^ + ln((v/^+l)2) 

S^ir) = -r-2v/f-ln((v^+l)2) (29) 
Next, having found the points Xi and X2, write the field at (A, r) as 

0(A, r) = - (01 (A + S^{r)) + 02(A + S2{r))) (30) 
r 

where the field on the initial surface is given in terms of two functions 

0(0, r) = i(0,(5,(r)) + 02(^2(r))) 
r 

= I (Mr) + Mr)) (31) 

Finally, in analogy with the fiat space case, write /i and /2 in terms of 0(r) and 7r(r) 
on the surface of quantization (A = 0). 



9 Black Hole: Thermometer Redux 

Let us now consider adiabatically switching on a thermometer, kept at fixed Schwarzschild 
r, and then switching it off. As in the case of the moving mirror, second order per- 
turbation says 

V{E,Eo) = ^\{E\M{Q)\Eo)\^ j dXdX 

^-i{E-Eo){\-\')^-{\-\of/2a 

e-^''-'°)'/^'^(0(A,r)0(A',r)) (32) 

Of course, now the points (A, r) and (A', r) in the correlation function (0(A, r)0(A', r)) 
are traced back, using nuU-geodesics, to the initial surface A = and the fields are 
appropriately re-expressed in terms of and tt on that surface. The only subtlety in 
this calculation is that for arbitrary r the interaction term gets an extra correction for 
the time dilation at point r, since the energy levels of the thermometer are defined 
in its rest frame. The result of this calculation is that the thermometer reads a 
temperature 

keT^ , ^ (33) 

SttM^I - 2M/r 

which agrees with Hawking's result. 
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10 Black Hole: Energy Flux 



The energy flux for the Schwarzschild black hole is calculated in the same way as for 
the moving mirror. Again we have to point-split the fields appearing in the energy 
momentum tensor and then take the limit of zero splitting. As before wc find that 
the flux Tx^r is finite and the total flux through a sphere at large r takes the expected 
limiting value 

plus, of course, transients which vanish for large A and terms which die faster than 



11 Back Reaction 

Our approach is unique in that we start at a finite time and calculate everything as 
a function for all values of A and r. Thus, in principle, one can discuss the problem 
of back reaction after the Hawking radiation has set in, but before any appreciable 
amount of the black hole's mass has been radiated. The reason there is a back reaction 
problem is because we calculated the energy-momentrum tensor for a Schwarzschild 
background and found a non- vanishing Txr and so we are in the situation that 

T^, ^ 0. (35) 

This, of course, is not consistent with the Einstein equations and so we do not have 
a self-consistent semi-classical problem. 

Given our computational procedure, however, for any point r, Txr is zero until 
radiation from the horizon has a chance to reach that point, at which time an observer 
begins to see the Hawking radiation. In principle we could feed our expression for 
T^iy back into the Einstein equations and attempt to find a self-consistent metric for 
which the computation of T^j, wouldn't change very much. Clearly this is difficult to 
do in general but we can ask what things look like inside of a sphere of radius r after 
the Hawking radiation has set in. If, in this region we adopt a metric of Schwarzschild 
form but change M to M{t) — Mq — Ft then the metric becomes 



+r^dn'^ (36) 
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which leads to an Einstein tensor of the form 





2F 





2F 





V 



where 






X{t) = 1 














4F^sin(6>)^ 



(37) 



2 (Mo -Ft) 



(38) 



which matches the outgoing flux computed for the static background. Since our 
computation of the flux at a point (t, r) only involves the computation of geodesies 
leaving the initial surface and arriving at this point. It is clear that for a large 
black hole these geodesies will not change very much for time intervals for which 
the Hawking radiation has set in but during which only a neglible fraction of the 
black hole mass has been radiated away. Therefore, it would seem that an iterative 
self-consistent solution should be possible. 



12 Entropy? 

While the two issues of entropy and the issue of decoupling of modes at r = are 
very interesting, there is no time to discuss them in this talk. A discussion of these 
issues will appear in a forthcoming paper. 



References 

[1] Kirill Melnikov and Marvin Weinstein, |hep-th/0109201| (2001), 

[2] S. W. Hawking, Commun. Math. Phys. 43, 199 (1975); J.B. Hartle and S.W. 
Hawking, Phys. Rev. D13, 2188 (1976). 

[3] W.G. Unruh, Phys. Rev. D14, 870 (1976). 

[4] T. Jacobson and D. Mattingly, Phys. Rev. D61, 024017 (2000). 

[5] M. Parikh and F. Wilczek, Phys. Rev. Lett. 85, 5042 (2000). 

[6] If the initial state is not the vacuum state but any finite energy state of the instan- 
taneous Hamiltonian, the finite energy will be radiated away in finite amount of 
time. The structure of radiation at large later times is independent of the initial 
quantum state of the massless field. 

[7] N.D. Birrell and P.C.W. Davies, Quantum Fields in the curved space, Cambridge 
University Press, 1984. 



11 



[8] B.S. DeWitt, in General Relativity, eds. S.W. Hawking and W. Israel, Cambridge 
University Press, 1979. 



12 



